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Measurement of high order current correlators. 
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a L.D. Landau Institute for Theoretical Physics RAS, 
Moscow, 119334, Kosygina str. 2 and 
b Theoretische Physik, ETH Zurich, CH-8093 Zurich, Switzerland 

The feasibility of measuring high-order current correlators by means of a linear 
detector is analyzed. Two different types of measurements are considered: mea- 
surement of fluctuation power spectrum and measurement of unequal-time current 



correlators at fixed points in time. In both cases, formally exact expressions in 
terms of Keldysh time-ordered electron current operators are derived for the de- 
tector output. An explicit time ordering is found for the current correlators under 
<— j , the expectation operator used in measurements of high order unequal-time current 



correlators. The situation when a detector measures current correlators at different 
points of a conductor is considered. 

PACS numbers: 03.65.Ta, 05.40.-a, 05.60.Gg, 73.23.-b 



I. INTRODUCTION 



Study of noise has recently been gaining popularity. This interest was quintessen- 
tially expressed by R. Landauer in the title of his paper The Noise Is the Signa^. 
It turned out that measurement of noise (electron current fluctuations) can be used 
to determine characteristics of electron transport that cannot be obtained by mea- 
suring only the average current or the conductance of the system. 

The phenomenon of shot noise was predicted by W. Schottky as early as in the 
beginning of the past century 2]. Schottky showed that the variance of the fluctu- 
ating current carried by statistically independent discrete charges is proportional to 
the product of the carrier charge with the average current. This fact can be used to 
measure the carrier charge in various systems. Originally, shot noise was measured 
in electron vacuum tubes [3J. Much later, measurements of shot noise were used to 

n n 

determine the quasiparticle charge in edge states for quantum Hall systems [4], [5[ 
and the double electron charge for Cooper pairs in hybrid NS systems [6j. Further- 
more, shot noise measurement can be used to study the electron correlations due to 
fermionic statistics and electron-electron interaction effects (e.g. see rewiev 

Considerable interest in study of current correlators was also motivated by 
prospects of preparation of bipartite entangled states in superconductors [8] and nor- 
mal electronic conductors [9], as well as of entanglement in quantum Hall effect 10]. 
The corresponding bipartite Bell inequality (which characterizes the degree of en- 
tanglement) is formulated in terms of second-order current correlators. Moreover, 

n fin 

measurement itself can be a source of entangled electron pairs [9j. In [U|, 112], a 



scheme was proposed for preparing an arbitrary n-electron entangled state, which 
requires measurement of nth-order current correlators. 

Furthermore, measurement of second and higher order current correlators can 
provide as much information as possible about the state of a conductor under con- 
straints imposed by quantum theory. This idea can be implemented in study of 
the full counting statistics of transmitted charge. Quantities of this kind have long 
since been analyzed in quantum optics in studies of the counting statistics of pho- 
tons emitted by various sources (e.g. see 13]). However, a similar idea regarding 

n n 

electron systems was proposed only in relatively recent studies \14L 1151 ] . where the 



distribution function was found for the charge transmitted through a quantum point 
contact over a finite time interval. 

The past decade has seen a rapid increase in the number of theoretical studies fo- 



cused on high-order current correlators, in particular in diffusive 16j,ll7(], chaotic 18] 



and interacting |19|, |20|, l2l(] e 



to the effects of temperature 



ectron systems. Considerable attention was also given 



22 ] and electromagnetic environment 23j on the third- 



order current correlator. A general approach to full counting statistics of current 
fluctuations, based on the path integral method, was recently developed in 



20 



24]. 



Despite numerous theoretical studies, only zero- or finite- frequency shot noise has 



y,y, 



26j). In 3, electron 



so far been amenable to direct measurement (see 
antibunching was studied in a fermionic Hanbury-Brown-Twiss-type experiment by 
measuring electron current cross correlations in a multiterminal conductor. Only in 
recent experimental studies, the third cumulant of electron current fluctuations was 
measured at low frequencies in a tunnel junction 28|, |29( and in a quantum dot in 
the Coulomb blockade regime 30]. 

These measurements of the third cumulant stimulated theoretical studies where 
various quantum detectors of current fluctuations were proposed. They can be di- 



vided in two types: (1) detectors where transitions between discrete energy levels 



are induced by current fluctuations 



31 



33 



34 . |35| |; and (2) threshold detec- 



tors where transition from a metastable state is caused by interaction with current 



fluctuations 



36 



32| • 



Since current operators taken at different times do not commute, practical mea- 
surement of current noise leads to an additional question about which particular 
correlation function is measured in an actual experiment. For second-order cur- 
rent correlators, the standard prescription is to use the symmetrized correlator 

391 ] the symmetrized correla- 



{I(t 2 )I(ti) + J(ti)J(t 2 )> [38|]. Indeed, as shown in 
tors at different times can be observed. Note, that at finite detector temperature, 
the anti-symmetrized correlator z(/(t 2 )-f (^i) — I(ti)I(t2)) also contributes to the de- 
tector output due to the back action of the detector on the measured current. An 
analysis of measurement of noise power spectrum by means of a resonant LC circuit 



coupled to the conductor was presented in 



40j (see also 3l|, |4jJ ) . In contrast to the 



measurement of the time resolved correlators, it was found that the ground-state 



(passive) detector can measure only the positive-frequency noise power spectrum 



where the current correlator is not time ordered. If the detector is in an excited 
state, then the negative-frequency power spectrum also contributes to the result. 
This behavior is explained by the fact that the ground state detector can only absorb 
energy, whereas an excited-state detector can transfer energy to the conductor as 
well. 

An analogous time-ordering problem for current operators under the expectation 
operator arises with regard to measurements of hi gh-o rder current correlators. In 
theana^es o f fu„ current stages presented inflQ, it was shown that the 
transmitted charge cumulants ((Q n (t))) with n > 3 depend on the current measure- 
ment method. In particular, when Q(t) is defined as the integral Q(t) = J * /(t')^' 



without time ordering, as in the naive approach developed in |14j, the resulting 
statistics corresponds to a fractional value of transmitted charge. However, an anal- 
ysis of a more realistic measurement scheme using an auxiliary spin- 1/2 system flo1 |. 
revealed that the charge counting statistics (far from the scattering point at x 3> x c , 
where for energy independent scattering x c is of the order of the Fermi wave length) 
is described by a binomial distribution corresponding to an integer value of the 
charge. In the latter scheme, the cumulants ((Q n {t))) are expressed in terms of 
Keldysh time-ordered electron current operators. This discrepancy between current 
statistics raised a question about which correlator can be measured in an actual 
experiment. 

In this paper, we theoretically analyze the feasibility of determination of high- 
order current correlators by measuring both fluctuation power spectrum and current 
correlators at fixed points in time. The present analysis does not address the dy- 
namics of electrons in the conductor and behavior of current fluctuations per se. We 
consider the joint evolution of the conductor and the detector and find expressions 
for the detector output in terms of correlation functions of current operators, which 
can be used to determine all measurable correlation functions of current fluctuations. 
Since the arbitrary correlation functions of general form are considered, the results 
obtained here apply to measurement of observables of any quantum system coupled 



to a detector. However, the analysis below is focused on the current measurement. 



Following 



39 



40| | we consider harmonic oscillator as a model of the detector coupled 
to the conductor at one or several points. In the general case, measurement of this 
kind corresponds to current measurement with a detector dynamics are described 
by a linear equation of motion. 

II. GENERAL SCHEMES FOR MEASURING CURRENT 

FLUCTUATIONS 

In this section, we describe two different schemes, which measure the power 
spectrum of current fluctuations at a finite frequency and current correlators at 
fixed times, respectively. We tentatively treat the detector as an arbitrary quantum 
system coupled to the conductor, assuming that measurement of its coordinate x 
provides information about the magnitude of the current carried by the conductor. 
In the sections that follow, the role of a detector is played by a resonant LC circuit 
or a single-electron ammeter. Since either device can be modeled by a harmonic 
oscillator, the problem technically reduces to analysis of a harmonic oscillator driven 



by an external field (e.g. see 43]), which provides formally exact expressions for the 
detector output. 

A. Continuous Quantum Measurement 

Suppose that the coupling between the detector and the conductor is adiabatically 
switched on at an instant to i n the past. At t > to, a measurement is performed on 
the detector observable x. The outcome of a quantum measurement on x can be 
described only statistically in terms of the distribution function of the value of x, 
or, equivalently, in terms of moments of the form (x n ), n G Z. 

This procedure is a typical example of continuous quantum measurement: the 
detector continuously interacts with the measured system, and the detector output 
is read out either at several points in time with limited accuracy or only once (see 



review m 



44] and references therein). The outcome is a time-averaged value of 



an observable, with an averaging kernel determined by the detector dynamics. In 



particular, if the detector is modeled by a harmonic oscillator, then a cumulant ((x 11 )) 
of lowest order in the coupling strength is proportional to the Fourier transform of an 
nth-order unequal-time current correlator or, equivalently, to the fluctuation power 
spectrum at the oscillator frequency (see below). 
The generating function of x defined as: 

x(X)=T rdst {p(t)exp(iXx(t))}, (1) 

where p(t) is the detector density matrix and the trace is taken with respect to 
detector degrees of freedom. Once x(A) is known, one can find (x n ) by differentiating 
X(A) with respect to A: 

<^) = HrHm(^ X (A)), (2) 

The generating function is determined by the detector density matrix p, which 
is obtained by tracing out the degrees of freedom of the measured system from the 
total density matrix D: p(t) = Tr sys {D(t)}. The density matrix D(t) is determined 
by the unitary evolution of the entire system starting from some point in the past. 
Suppose that the joint density matrix of the detector and the conductor at the 
starting point t is a direct product of the form D(t ) = pi n Ri n , where the density 
matrix Ri n represents the initial state of the conductor. We introduce a Hamiltonian 
H sys to describe free dynamics of the conductor and a Hamiltonian H de t = Ho + Hmt, 
to describe detector dynamics, where H is the free detector Hamiltonian and H ult 
is the coupling Hamiltonian. In the interaction representation with respect to H sys , 
the total density matrix at t > t has the form: 

D[t) = S(t , t) Pin R in (t)&(t , t), (3) 

where the operator S(t ,t) describes the evolution of the detector under the Hamil- 
tonian -f^dct : 

S(t ,t) = Texp(-^J^H dct (t')dt'), (4) 

with T is the forward time ordering operator. 

We change to the x representation for the total density matrix and the system 
evolution operator: 

D(x,y,t) = (x\D(t)\y), 



S(x ,t ;x,t) = (x\S(t ,t)\x ), 

which remains operator acting on the state of the conductor. Now, the characteristic 
function x(A) can be represented as: 

X(A) = Tr sys ^R in (t) J dxe lXx p in (xo,yo)S ] (yo,t ;x,t)S(x ,t ;x,t)Y (5) 

For A = 0, the expression on the right-hand side here is the influence func- 
tional [431 ] . describing the back-action of the detector on conductor dynamics. 



B. Subsequent Quantum Measurements 

Now, consider another feasible measurement scheme, where the detector state 
is measured at subsequent time moments t\ < t 2 < ■■ < t n . According to von 



Neumanns projection postulate 45(, an ideal quantum measurement performed on 
x at any £j projects the detector state at ti onto an eigenstate of the operator 
x. This measurement is strong in the sense that it strongly changes the detector 
wavefunction at point ti and therefore affects its further evolution and measurement 
outcomes at t > ti. 

First, consider measurements performed on x at two subsequent moments t\ and 
t 2 . The resulting correlation function (x(ti)x(t 2 )) can be expressed in terms of the 
joint probability P(xi, t\\ x 2 , t 2 ) that the values of x measured at ti and t 2 are x\ 
and x 2 , respectively: 



{x{ti)x{t 2 )) = J dxidx 2 xix 2 P(xi,ti,x 2 ,t 2 ). (6) 
In turn, the joint probability P(x%, t\\ x 2 , t 2 ) is given by the product: 

P(xi, ti] X 2 , t 2 ) = P(x 2 , t 2 \x!, tx)P(xi, h), 

where P(x 2 , t 2 \x\, ti) is probability that the value of x is x 2 at t 2 conditioned on the 
outcome x\ of the measurement at moment ti, and P{x\,t\) is the probability to 
measure the outcome x\ at moment t\. 

The last two probabilities can be determined from principles of quantum mechan- 
ics. If D(t) is the joint density matrix of the detector and the conductor, then 

P(xi,t 1 ) = Tr 8y8 {(xi|D(t 1 )|x 1 >}. 



The conditional probability can be found by invoking von Neumanns projection 
postulate: at time moment t 1 , the total density matrix undergoes an instantaneous 
change: 

h(+-0\ ^ f>(++o\ _ ki)(sil-P(*r°)ki)(sil m 
Tr s ^{{x x \D{t x )|aji)> 
where the denominator is introduced to normalize the reduced density matrix. At 

t > ti, the reduced matrix is described by a unitary evolution operator: 

D(t 2 ) = S(t 1 ,t 2 )D(t+°)S\t 1 ,t 2 ). 
Thus, the conditional probability P(x2, t 2 \xi, h) is expressed as follows: 

P(x 2 ,t 2 \x 1 ,t 1 ) = Tr sys \S(x 1 ,t 1 ;x 2 ,t 2 ) — E^iEll^L & (x l ,t 1 ; x 2 ,t 2 )\, (8) 

Tr sys {/;(a;i,a;i,ti)} 

where D(x, y, t) = (x\D(t) \y) is the total density matrix in the basis of eigenfunctions 
of the observable x. 

Taking the product of the conditional probability on P(xi,ti), we find that the 
normalization factor cancels out. The result is 

P(x ± , h] x 2 , t 2 ) = TrsysiSix!,^ x 2 , t 2 )D(x 1 ,x 1 ,t 1 )S\x 1 ,t 1 ] x 2 , t 2 )}. (9) 

Since the correlation function (x(ti)x(t 2 )) is calculated by averaging over all pos- 
sible outcomes at t\ and t 2 , the measurement performed at t x can be interpreted as 
the instantaneous diagonalization of D(x,y,t) in the basis of the measured observ- 
able, 

D(x,y,t^)^D(x,x,tt ), (10) 

which is an equivalent formulation of von Neumanns projection postulate. 

In the general case of subsequent measurements, the correlation function of the 
detector outputs at times tj can be calculated as: 

/n 
YlxidxiPdxuU}), (11) 
i=l 

where P({xi, ti}) is the joint probability that the detector coordinates at ti, t n are 
x±,...,x n respectively. For a single measurement of the detector coordinate x, the 
probability P(xi,ti) is expressed in terms of the diagonal elements of the detector 
density matrix at time t\: 

P(x 1 ,t 1 ) = Tr sys {D(x 1 ,x 1 ,t 1 )}. 



By analogy, the probability P({xj,tj}) for subsequent measurements can be ex- 
pressed in terms of diagonal elements of the density matrix 



D({xi,Xi,U}) = J dxodyoS(x n -i,t n -i;x n ,t n ) ... S(x ,t ;x 1 ,t 1 ) pi n (x ,yo)Ri n 

x &(y 0l to;x ll t 1 )...Sfl(x nr .i,t n - 1 ;x n ,t n ), (12) 

where the values of x at t — t i have to be equal to each other under forward 
and backward time evolution. Thus, the detector density matrix at the times of 
measurements, t = ti, is diagonal in the basis of eigenstates of x. Therefore, the 
required probability is 

P({ Xi , U}) = Tr sys {b({ Xi , Xi , ti})}. (13) 

Note that the correlation function defined in this manner must be real-valued, 
because each measurement outcome is a real number by virtue of von Neumanns 
projection postulate. This correlation function should be distinguished from the 
correlation function (x(ti)...x(t n )) of coordinate operators, which may not be real- 
valued in the general case. 

In this approach, each measurement projects the state of the system onto an 
eigenstate of the measured observable, as in von Neumanns ideal measurement. 
However, it is not the only type of feasible time resolved measurements. For example, 
the theory of photon detection in quantum optics assumes that the photodetector 
is reset to the ground state after each count and then switched back to the standby 

n 

mode [13J. In distinction to (JTj), this measurement of x(ti) = Xi formally corresponds 
to an instantaneous transformation of the density matrix bringing the system into 
the next state, 

D(tr) - D(t+) =Pi® foWpl-E*} ? (M) 

Tr sys {(xi|D(tj )\Xi)} 

which corresponds to the direct product of the detector density matrix pi with the 
conductor one. The density matrix pi represents the detector state in standby mode 
at time U. The corresponding correlation function of x is expressed as: 

(x(t 1 )...x(t n )) = Tr det {x S{t n -i,t n )p n - 1 x ... x Tr det {x S(t 1 ,t 2 )pi 

x Tr det {x S(t , tJptoRin&ito, ti)}^*!, t 2 )} x ... x 5t(t„_ a , *„)}• (15) 



In a sense, this type of measurement is more natural than von Neumanns ideal 
measurement. Indeed, the procedure described here assumes that the detector itself 
is a quantum-mechanical system whose state must also be measured. The state 
of the detector can be measured by coupling it to some macroscopic device. The 
coupling must be sufficiently strong for the detector state to be determined over 
a time interval much shorter than the interval between subsequent measurements. 
The resulting detector state may be different from |a^}, for example, in the ground 
state. This occurs when a detector excited to a metastable state by interaction with 
the measured system loses energy to the macroscopic device. 



III. MEASUREMENT OF NOISE POWER SPECTRUM 

Consider a measurement of the power spectrum of current fluctuations using 
a resonant LC circuit inductively coupled to a quantum conductor as a detector. 
Suppose that the observable to be measured is the electric charge q{ti) stored in 
the capacitor at a point in time t\. As a model of the LC circuit, consider a weakly 
damped harmonic oscillator driven by an external force J(t) = J2i oiiIi(t). The total 
Hamiltonian is 

#d* = -!" + ^+JE«i£(*), ( 16 ) 

where L and C denote the inductance and capacitance of the circuit, the operators 
j and Jj represent the currents through the LC circuit and through the conductor 
at coordinate xf, and a« denotes the mutual inductance between the LC circuit and 
the conductor in the neighborhood of X{. The operators q and j can be expressed 
in terms of bosonic creation and annihilation operators: 

^(iac) (a +a)) 3 '(21) (a a)) (17) 

where f2 = 1/V LC is the resonant frequency of the LC circuit. 

For the harmonic oscillator with Hamiltonian ([161) . the evolution operator in the 
q representation can be found in explicit form: 

S(q h q 2 ,T) = S ( qi ,q 2 ,T)TexJjr^^l^ j J(t) cos n(f a -t)dt (18) 



_ T / ^C*) 008 ^ ~ tl " )dt+ j2 J dt j dsJ(t)J(s) cosQ{t 2 -t)cosfl{s -*i)} J, 
*i *i ti ' 

where T = t 2 — t\, and the transition amplitude for a free harmonic oscillator is 

expressed as 

Assuming that the LC circuit prior to measurement (at t — > — oo) is in equilib- 
rium at temperature and and performing the integral in (j3j) over the initial and 
final detector coordinates, we find that x(A) is the product of the equilibrium char- 
acteristic function Xo(A) with the characteristic function of excess capacitor-charge 
fluctuations: 



X(A) = Xo(A) (T ± exp 



F[J + J- 



l — ( 1 (G(ti - t)JJt) + G*(ti - t)J_(t)) dt 
n Jt 

'(20) 

where the influence functional of the detector F[J + , J_] is expressed as 

F[J + , J_] = exp ~ f 1 dt f ds (J + (t) - J_(t))(T(t - s)J + {s) - T*{t - s)J_(s))] , 

(21) 

the time-ordering operator T± is introduced to order the current operators J±(t) 
forward or backwards in time, and averaging is performed over the detector states. 
The response functions in (I20p are defined as follows: 

G(t) = (A)((jv + l)e~ int - Ne int ) (22) 
T(t) = (—) {(N + l)e~ im + Ne int ), (23) 



2L 

where N = l/(e m/e - 1) is the bosonic occupation number for the oscillator mode. 
The equilibrium characteristic function XoW corresponds to Gaussian charge fluc- 
tuations with variance Sq 2 = (h/LQ)(N + |): 

Xo (A) = exp(-~AV). (24) 

Expressions (1201) and ( 12~TT) can be rewritten in more compact form in terms of the 
Keldysh Greens functions 

G K (t-t') = -i(T K q(t)j(t')), 



r K (t-t') = (T K j(t)j(t')), 

where the operator denotes time-ordering along the Keldysh contour consisting 

of forward and backward branches (ioj^i) an d (^i,^o) : 

"zA 
h Jk 



xW = Xo(A)(7ft:exp 
1 



x cxp 



2K 1 



K 



G K (ti -t)J(t)dt 
T K (t- s)J(t)J(s)dtds 



(25) 



This expression provides a formally exact description of capacitor charge statistics 
for the LC circuit. The characteristic function x(A) can hardly be calculated in 
explicit form, but expression (1251) makes it possible to develop a perturbation theory 
in the parameter oc/L <C 1 by representing the influence functional F[J + , J_] as a 
Taylor series expansion. 

The analysis that follows focuses on the limit case when the back-action of the 
detector on the measured system, described by the influence functional given by ( 12~TT) . 
is negligible. Setting F[J + , J_\ = 1, we obtain: 



x(A) = Xo(A) (T^exp 



*4 I G K (ti-t)J{t)dt 



(26) 



First, we examine the autocovariance of the capacitor charge fluctuations. Us- 
ing ([2]), we obtain: 



(27) 



where the first and second terms correspond to equilibrium and excess fluctuations, 
respectively. 

It should be noted here that perturbation series expansions contain secular 
terms [46]. Indeed, even though the exact expression for x(A) given by eq. ( 1251 ) 
formally yields a finite result for any A, each term in the perturbation series gen- 
erally diverges as t\ — ► oo. This is obviously true even for second-order quantities: 
in the stationary case, the current-current correlator in (j2*7|) is determined only by 
a time difference, and the resulting integral diverges as one of the remaining time 
variables goes to infinity. 

These divergences can be eliminated by using weakly damped Greens functions. 
In physical terms, this means that a real detector not only interacts with the con- 
ductor, but also loses energy to the environment. 



However, the introduction of an external heat bath is not a physical necessity. 
Indeed, its role can be played by the detector itself. The process can qualitatively 
be described as follows. An interaction with the conductor brings the detector into 
a steady state that can be strongly overheated as compared to its initial state. Any 
further excited state of the detector has a finite lifetime, because the detector will 
tend to lose energy to the conductor. 

The degree of detector overheating can be estimated perturbatively. For the 
detector starting from the state with occupation number N, we use perturbation 
theory to calculate the probabilities of excitation to the state N+l and de-excitation 
to the state N — 1 over a time interval At ^> Q^ 1 : 

P N+ i(At) = ^^+l)S (2) (fi), (28) 
Pn-i(M) = ^^NSW(-n). (29) 

In perturbation theory, these probabilities are low. The power spectrum of current 
fluctuations is defined as 

4 2) (u,) = J (Uh^ih)) e-^-^dih - t 2 ). (30) 

Using the balance condition Pn-i > Pn+i, we estimate the steady-state occupa- 
tion number as 

This quantity may be large even in the case of weak coupling between the detector 
and the conductor, which becomes obvious as the excitation energy tends to zero 
(fi -> 0). 

However, we assume below that the detector is coupled to the environment, in- 
troducing the factor e _r? '*' into integrals in all expressions analogous to (|2Tj) . To 
estimate the leading-order contribution to the excess capacitor-charge variance (q 2 ), 
we express the average over detector states in (|27p in terms of integrals of the noise 
power spectrum: 



In the weak-damping limit (77 <C Q), we have r]/(x 2 + rj 2 ) — > 7r5(x), and the result 
is expressed in terms of the positive- and negative- frequency power spectra, 5^' (fi) 
and S^\-Q): 

<? 2 > = E SI 1 W(n) + N{s§\n) - sf{-n))}. (33) 

* J I — I 
Expression (1331) extends the results presented in [401] to the case when the detector 

is coupled to the conductor at several points. The result for current measurement at 



two different points was presented in 



471 ] and calculated later perturbatively in 



It is clear that negative-frequency current fluctuations contribute to this result 
only at finite temperatures (N > 1), whereas positive-frequency fluctuations con- 
tribute even at zero detector temperature. As noted in 40, |4l| the reason is that 
energy transfer at zero temperature is possible only from the conductor to the de- 
tector, and the intensities of such processes are proportional to Sj^(f2). At a finite 
temperature, the detector can be in an excited state from the start and transfer 

(2) 

energy to the conductor at a rate proportional to SL (— fi). 

The proposed method can be used to measure any electron current correlation 
function. However, fluctuations of capacitor charge are dominated by fourth-order 
correlators. Indeed, it can easily be shown that the lowest order contribution to the 
mean stored charge, 

<»4?£®V<>- <*) 
is small since rj/Q <C 1. An analogous result holds true for any contribution to 
charge fluctuations containing odd-order correlators. 

Now, consider measurement of fourth-order correlators. Retaining only the lowest 
order nonvanishing contribution to the fourth-order cumulant ((? 4 )) = (q 4 ) — 3(q 2 ) 2 , 
and using characteristic function (1261) . we obtain: 

«g 4 » = ^(T K (J K G K (t 1 -t)J(t)dt)^-3(q 2 ) 2 . (35) 

For simplicity, suppose that iV = at the starting point. Then, the dominant 
contribution to capacitor charge fluctuations corresponds to the following time or- 
dered average of current operators: 
fi °r 

((<? 4 » = j^- A j ds l ds 2 dt 1 dt 2 e in ^ +t ^-^ (T^{J(s 1 )J(s 2 )}T + {J(t 1 )J(t 2 )}) . 



— oo 



(36) 



The fourth-order unequal-time current correlator contained in this expression can 
be represented as the sum of an irreducible correlator and a number of reducible 
ones: 

(W3/4) = {{hkhh)) + (/i/ 2 >(/ 3 /4> + (ileitis) + {hh){hh) + 

where all odd-order correlators are dropped. Since steady-state current correlators 
depend only on time differences, the spectral density of the fourth-order irreducible 
correlator can be represented as 

S\f kl {u U LO^) = J nrf(t i -t m )e-^^-^ 1 )((/ J (t 1 )/ i (t 2 )4(t3)/K^))), (37) 

8=1 

whereas the reducible ones can be expressed in terms of the noise spectral den- 

( 2) _____ 

sity S!-j (u) (cf. eq. (13111) ). Introducing a damping factor and performing the time 



integration, we obtain: 

((<? 4 » = 24£ 



aiOtjakai 



X 



X • 



+ 



ijkl 

7 W 



du 1 du 2 du 3 S ijkl (ui,uj 2 ,uj 3 ) 

(2tt)3 (ui-n- iri)((u 2 - 2Q) 2 + 4r] 2 )(u 3 -VL + irf) 

dujiduj 2 S\f{u x )S^{u 2 ) ( 1 



(2vr) 2 (ux +uj 2 -2n) 2 + 4^ 2 Ulu 2 - tl) 2 + V 



2 



+ 7 Fi "T7 Fi ^H-S^ 2 ) 2 . (38) 

Finally, we note that the output of a measurement of the power spectrum of 
current fluctuations performed with the same LC circuit at several points contains 
not only cross-covariances (SV, (±f2) for i 7^ j), but also auto covariances (Su(±Q) at 
the same point). To measure only cross-covariances, one must use two independent 
LC circuits with equal resonant frequencies coupled to the conductor at only one 
point. The results presented above can be extended to the case of several detectors. 

In particular, for two LC circuits coupled to the conductor at points x and y, 
respectively, it can be shown that the charge correlator (q x (ti)q y (ti)) calculated 
to second order of perturbation theory is proportional to the second-order cross- 
covariance of current fluctuations: 

M) = j0 ^{(n + i)(s£?(n) + - N(s^i-n) + s$(-n))}. (39) 



IV. MEASUREMENT OF UNEQUAL-TIME CURRENT 

FLUCTUATIONS 



Consider an experiment where current correlators are to be measured at certain 
points in time. In a thought experiment of this kind, we can use single-electron 
ammeter, whose readout is proportional to the exact instantaneous value of current. 
As a classical model of the ammeter, we can consider a damped harmonic oscillator 
driven by an external force proportional to current: 

x{t) + rjx(t) + tt 2 x(t) = a I{t), (40) 

where x(t) is the ammeter readout, Q is the resonant frequency, rj is the damp- 
ing factor, and a is the coupling constant. The solution to this equation can be 
represented in terms of the response function x(t) = J x(t — t')I(t')dt': 

f du . . ,■ t f duj -a e~ luJt , . 

J 2tt J 2tt ui — & 2 — | — iTjuj 

Since the oscillator motion is described by a linear equation, we can use the classi- 
cal dynamical susceptibility xi 00 ) °f the damped oscillator to determine temperature 
Greens functions by invoking the fluctuation dissipation theorem 49[]. Accordingly, 
we first express the correlation function (x(ti)...x(t n )) of detector outputs in terms 
of Greens functions for the free oscillator and then replace these functions with those 
for the damped oscillator. 

Thus, we assume that the detector output is the ammeter readout x(t), and the 
ammeter is described by the Hamiltonian: 

H det = ^- + — amxl(t), (42) 

2m 2 

where the mass m and the momentum operator ft play the roles of moment of inertia 
and angular momentum, respectively. Unlike LC circuit Hamiltonian ( fl6l) it is linear 
in the coordinate rather than momentum of the oscillator. 

In the coordinate representation, the amplitude of the transition from the state 
(xi,t]) to the state (x 2 ,t 2 ) driven by the external current I(t) is expressed as fol- 



lows 



43|: 



S(xi,t 1 ;x 2 ,t 2 ) = S (x 1 ,x 2 ,T) Texpl in ^ T j J I(t)smSl(t 2 - t)dt (43) 



t2 2 *2 t 

+ ~~TT I ^ sin n(yt ~ tl ^ dt ~W I dt I ds sin ~~ ^ sin n ( s ~ tlj 



ti ti ti 

where S (xi, x 2 ,T) is the transition amplitude for a free oscillator given by (TH31) . 

As shown above (see eq. (1121) ) , in the case of an ideal measurement of the coordi- 
nate x(ti) at each point in time, the correlation function (x(ti)...x(t n )) is calculated 
as the average over the intermediate coordinates of evolution operators of the form: 

where Sij = S(xi, U\ Xj,tj). However, it is clear from expression (l4"3"j) for the evolution 
operator that such an average is not well defined. Indeed, as a function of Xi this 
product of evolution operators contains the factor: 

T±exp — / J + * -/_ *))_ -A 

+ r i (/,w-f-W) f"^ 1 "* ) , 4)- (44) 

sinS2(ti+i - ti) ) j 

where 7± denotes the time-ordering operators introduced to order the current op- 
erators I+(t) and forward and backwards in time, respectively. The current 
operator under the time-ordering symbol can be treated as ordinary functions of 
time. This demonstrates that a weighted average of this form of the coordinate Xj 
is a divergent function of I±{t). 

One possible explanation of this divergency is that the oscillator with exactly 
known coordinate Xi has infinite energy. Accordingly, the probability that the oscil- 
lator will have a certain coordinate Xi + \ at a subsequent point in time t is indepen- 
dent of Xi + \ 



[e.g. see 



2nh sin Qt 

(431). 

This unphysical result can be avoided by several means. First, we can assume 
that the coordinate is measured with limited accuracy \x — Xi\ ~ A. Then, the 
postmeasurement state of the oscillator is described by a density matrix localized 
within the measurement error, such as p Xi (x,y) ~ exp(— [(x — Xi) 2 + (y — Xj) 2 ]/4A 2 ). 
Second, as noted above, the measurement maps the detector density matrix to a 
certain standby state. 



Hereinafter, we consider the latter scenario and assume that the measurement on 
x(ti) brings the detector into an equilibrium state at temperature 0. Then, the cor- 
relation function of detector outputs is given by formal expression ( fl~5l) . Integrating 
the evolution operators over the initial and final coordinates, we obtain: 

{x{h)...x{t n )) = lT K F K [i]f{Q J C K {U-t)i{t)dt)\ (45) 

I — 1 

where the Keldysh Greens function is defined as: Cxit — s) = (TK{x(t)x(s)}); 
K\ = [ti + iO — > tj —>■ ti — iO] - is the Keldysh contour. The influence functional 

f k [i\. 



F K [I] = exp 



2h 2 



C K {t- s)I{t)i{s)dtds 



K 



(46) 



is defined on the Keldysh contour K c _ 

The Keldysh Greens function Cx{t) can be expressed in terms of the causal 
Greens function C(t > 0) = (x(t)x(0)), which can in turn be expressed in terms of 
the dynamical susceptibility given by xi^) pdj) - by invoking the fluctuationdissipa- 
tion theorem: 

<*> = I^e- W 

In what follows, Greens function C (t) is conveniently represented as the complex 
variable: C(t) = S(t) + iA(t), with real and imaginary parts corresponding to 
the symmetrized and antisymmetrized coordinate autocorrelation functions: S(t) = 
~(x(t)x(0) + x(0)x(t)), and A(t) = ±(x(t)x(0)-x(0)x(t)). Using the relation 1/(1 - 
e~ x ) = | + | coth(x/2), we obtain: 



S(t) = l^J d u cosM) lm X {uj) coth(^) , (48) 
A(t) = / du sin(cjt) Imx(a;). (49) 

27T J 



A. Measurement of Second-Order Current Correlators 

Here, we derive an expression for the correlation function (x(t\)x{t2)) in a per- 
turbation theory in the coupling constant a < 1 thus neglecting the back-action of 
the detector on conductor dynamics (i.e., by setting Fk[I] — 1). The perturbative 



calculation of {x(ti)x(t 2 )) was first carried out in Ref. [39(] . Rewriting Eq. (T4"5"]) in 
explicitly time-ordered form, we find 

2 * 2 ' J 

(s(fi)a:(f 2 )) = -|r/ < {A(f a - ti)</(f 2 )fyi)> S 

-L4(t 2 - 4)5^ - t' 1 )(/(t 2 )/(t' 1 )) A }, (50) 



where (I(ti)i(t 2 )) s and (/(ii)/^)}" 4 are the symmetrized and antisymmetrized cur- 
rent correlators, 

(m)Hh)) s = \m)i(t 2 ) + i(f a )i(fi)>, (51) 

<J(t 1 )/(^)) A = ^(/(ti)/(t 2 ) - I(*2)/(*i)>. (52) 

Since we focus on the measurement of current correlators at fixed times, t\ and t 2 , 
the functions A(t) and S'(i) in Eq. (1501) must decay over time intervals much shorter 
than the time scale r e , of the evolution of the electron system. Otherwise, the 
correlator (x(ti)x(t 2 )) would be proportional to time-averaged current correlators. 
Calculating the integral over frequency in Eq. ( f49i) . we find that A(t) is a decaying 
oscillating function independent of temperature: 

A(t) = exp(-^-) sm{u v t), (53) 



where u v = yf2 2 — rf /A - is the renormalized oscillator frequency. 

The symmetrized correlation function S(t) is the sum of two functions, S(t) = 
Si(t) + S 2 (t), having different long-time behavior. The decay of S\(t) is completely 
determined by the oscillator parameters: 



s 1 (t) = -A. e -w2| e ^t coth + Cc x (54) 

whereas the decay of the other function depends on the oscillator temperature: 



aw = -^5 M + ^-^ - (55) 

where z/ n = 2imQ/h is a Matsubara frequency. 

Since the oscillator frequency is renormalized with Q held constant, the functions 
A(t) and Si(t) cannot decay faster than e~ nt (at 77 = 2fi) as i — > 00. In the 



strong-damping limit (77 ^> Q), the damping factor is smaller, and both A(t) and 
Si(t) decay approximately as ~ e ~ n ( n /v)t^ Therefore, measurements of unequal-time 
correlators should be performed with oscillator parameters chosen so that 77 ~ 2Q ^> 
max{r e _1 , \t 2 — tiT 1 }- 

Since the decay of S 2 (t) depends on the detector temperature, it is necessary 
that v\ ^> r" 1 , which is equivalent to the requirement that G 3> max{hr~ 1 , h\t 2 — 
til -1 }. For a ballistic conductor, the time scale of current fluctuations is r e ~ 
H/eV, where eV is the voltage across the conductor; therefore, the condition ;§> 
hr~ l corresponds to 6 > eV. This means that weakly nonequilibrium current 
fluctuations will be measured when the detector and the conductor have comparable 
temperatures. 

In this approximation, the current operators in (|5"U1) can be factored out from the 
integral, and we can set: 

£ 2 A(t 2 -t)I(t)dtn-^I(t 2 ). (56) 
£ S(t 2 -t)i(t)dta^I(t 2 ). (57) 

As a result, we find that the quantity measured at high temperatures is the following 
combination of the symmetrized and antisymmetrized current correlators: 

(sfaMfc)) = ^{<^)/(ti)> 5 + K^§) (i(t2)i(h)) A }. (58) 

At 77 ~ 2Q the relative contribution of either correlator depends on the ratio 40/M7, 
which can have any value in this regime. 

The high temperature result ( 1581) in the limit lu v <C rj was obtained in Ref. 39]. 
Another interesting experiment is measurement of quantum current fluctuations in 
the low-temperature regime (0 <C hr^ 1 ). In this regime, measurement of highly 
nonequilibrium excess fluctuations becomes particularly feasible. In contrast to the 
high-temperature regime, the long-time decay of S(t) follows a power low, 

S( t >n-) = -(^)l. (59) 

Again, requiring that 77 ~ Q ^> max{r e _1 , \t 2 — we can assume, 

j s{t 2 -t)i(t)dt^i(t 2 ) j dts{t)n±[Q + - r _ t) )Hh), (60) 



where in the last equation we have used the relation 

S{t)dt = r]e/Q\ (61) 



which holds at any temperature. The resulting expression for (x(ti)x(t2)) at low 
temperature of the detector takes the form, 

(x(h)x(t 2 )) = ^{(Ht2)Hh)) S + i^(Q + ^_ t J (/(t 2 )/(*i)) A } (62) 

Comparing the last expression with the high-temperature result (I58p we conclude 
that in the low-temperature limit the contribution of the antisymmetrized current 
correlator to the (x(t\)x(t2)) is proportional to the effective temperature, 

e' = e + ^, (63) 

where At is the time separation between subsequent measurement times. 



B. Measurement of Higher Order Correlators 

Perturbative treatment of Eq. (I45p makes it possible to express a coordinate 
correlation function of any order in terms of current correlators. First, consider 
measurement of third-order correlators. Time-ordering the current operators along 
the Keldysh contour, we obtain 

(xitJxfoMta)) = -( — ) / dtU dt' 2 dt[x 

v a ' Jt 2 h x Jt 

x{A(t 3 - t' 3 )A(t 2 - t' 2 )A( tl -t[) {I(foI&)i$ L )) 8 

-iA(t 3 - QA(t 2 -t' 2 )s( tl -t[) (iy a )i&)i&)) 8M 
-iA(t 3 -t' 3 )s(h- t' 2 )A( tl -t[) (iy a )i&)i&)) SAa 

+A(t 3 - QS(t 2 - t'JSih-t'J (i(t' 3 )i(t' 2 )i(t[)} A }. (64) 

It is clear that four different third-order current correlators are measured in this 
case. Denoting by {A, B} and [A, B] the anticommutator and commutator of opera- 
tors, respectively, we represent two of the four correlators as averages of symmetrized 
and antisymmetrized combinations of current operators taken at different times, 

<Wi> S = i({{/ 3 ,/ 2 },A}), (65) 
(Wi) A = ^([[/ 3 ,/ 2 ],/i]), (66) 



and the other two as averages of mixed combinations of symmetrized and antisym- 
metrized current operators, 

(Wi) 5Al = ^<[{/ 3 ,/ 2 },A]). (67) 
(Wi) 5A2 = ^({[/ 3 ,/ 2 ],A}). (68) 

Note that these four Hermitian combinations of current operators are not the only 
possible ones. For example, following a standard prescription, we might assume that 
the completely symmetrized combination is measurable: 

(A/2/3) 5 ' = l{h{hJ 3 }+Ukh}+UhJi}). (69) 

However, the analysis above shows that this is not true. The formal reason is that 
the current operators used to describe detector dynamics are Keldysh time-ordered. 
Since ti < t 2 < t 3 , the combinations (I(t 2 )I(ti)I(t 3 )) and (i(t 3 )i(ti)I(t 2 )) , which 
are contained in (1691) and not contained in (IBl)]) . can not appear as a result of such 
time-ordering. Indeed, the operator I{ti), to be measured first will be on the right 
or left of other current operators, depending on whether time- or antitime-ordering 
is performed, respectively. In other words, the absence of {I(t2)I(ti)I(t 3 )) and 
(/(t 3 )/(ti)/(t 2 )) in correlator fl65|) is dictated by causality: measurements performed 
at later points in time t 2 and £3 cannot affect the outcome of the earlier measurement 
at £1. 

This argumentation suggests that a linear detector can be used to measure only 
correlators of the form: 

(z(ti)-z(*„)) ~ ({{/(*„), (70) 

where t\ < ... < t n , as well as the remaining 2 n ~ 1 — 1 correlators obtained by 
replacing any number of anticommutators with commutators. 

Note that, depending on the measurement regime, the contributions of some of 
these 2 n_1 correlators to (x(ti)...x(t n )) may be equal. In particular, consider the 
measurement of a third-order correlator, 

Wj 1 ) I (j 2 )x(t 3 )) = ^{(/>3)/(i2)/(tl)> S + <^)(/(«3)/(i2)/(t,)> S ' 11 



where 6'i = 6 H — ,, - -. , and 6' 9 = 9 H — ,, a , , are two effective temperatures. In 

7T(ti— to) Z 7T(t2— tl) r 

the high-temperature regime one has Q[ ~ 0' 2 ~ G and hence, 

(x(t 1 )x(t 2 )x(t 3 )) = ^{(/(t 3 )/(t 2 )/(ti)) S 

+i^§) (/(t 3 )/(t 2 )/(ti)) SA3 + (H) 2 a(t 3 )/(t 2 )/(ti)) A }, (72) 

where the correlator {I^I 2 Ii) SAz is a combination of (i 3 i 2 Ii) SAl and (/a^/i) 5 ^ 2 
having the form: 

(Wi> 5A3 = ^ (<Wi> - (Ws>). (73) 



V. CONCLUSIONS 



We have theoretically analyzed measurements of higher-order current correlators 
with a quantum detector coupled to the conductor. A damped harmonic oscillator 
with coordinate x(t) measured directly is considered as a model of a detector whose 
classical dynamics are described by an equation of motion linear in x. 

Two essentially different types of measurement have been considered: measure- 
ment of the power spectrum of current fluctuations and measurement of unequal- 
time current correlators at several points in time. 

In the former measurement, a weakly damped resonant LC circuit inductively 
coupled to the conductor is used as a detector. In measurements of this type, 
the detector interacts with the conductor over a considerable time interval, and 
the measured observable q is the charge stored in the capacitor of the LC circuit. 
Measurement of moments (q n ) of the stored charge provides information about nth- 
order power spectrum of current fluctuations. A complete statistical description of q 
is given by Eq. ( 1251) in terms of integrals of Keldysh time-ordered current operators. 

This equation is used to develop a perturbation theory in the coupling strength. 
By neglecting the backaction of the detector on the measured system, expression 
( 13*31) is obtained for the second-order irreducible charge correlator ((q 2 )), in the case 
when the detector is coupled to the conductor at several points. At zero detector 
temperature, {(q 2 )) is proportional to the positive frequency Fourier transform of 
the second order current correlator symmetrized with respect to coordinates. At a 



finite detector temperature, negative-frequency current correlator symmetrized with 
respect to coordinates also contributes to the stored charge. 

Measurements of unequal-time current correlators are analyzed in the strong- 
damping limit, when the damping factor comparable to the oscillator frequency. 
The measured quantity is the correlation function (x(ti)...x(t n )) of the oscillator 
coordinate at subsequent points in time. The oscillator is assumed to relax to 
equilibrium after each measurement. This type of measurement is different from 
von Neumanns projective measurements at subsequent time moments. It is found 
that the correlation function (x(ti)...x(t n )) may be difficult to determine because an 
electron in a sharply localized state has infinite energy. 

In the general case, the outcome of time resolved measurements of the detector 
state is described by Eq. (j4"5|) . Perturbation theory is used to derive expressions (J5U1) 
and (1641) for second- and third-order coordinate correlators, respectively. When 
second-order correlators are measured, the detector output combines contributions 
of current correlators symmetrized and antisymmetrized with respect to time. In 
the high-temperature regime (0 3> eV) the contribution to {x(ti)x(t2}) due to 
the antisymmetrized current correlator is proportional to the detector temperature 
6, see Eq. (1551) . At low temperatures, (9 <C eV) one has the same result for 
(x(ti)x(t2)) as in the high-temperature regime with an effective detector temperature 
6' = Q+h/(irAt), where At is a time separation between subsequent measurements, 
see Eq. (1S2|). 

When higher order coordinate correlators are measured, a larger number of differ- 
ently time-ordered current correlators contribute to the measurement outcome. The 
current operators in such correlators are Keldysh time-ordered in accordance with 
the causality requirement: a future current measurement cannot affect the outcome 
of a past measurement. 
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